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The phase transition of the even offspringed branching and annihilating random walk is studied 
by iV-cluster mean-field approximations on one-dimensional lattices. By allowing to reach zero 
branching rate a phase transition can be seen for any TV < 12. The coherent anomaly extrapolations 
applied for the series of approximations results in vx = 1-85(3) and j3 = 0.96(2). 



I. INTRODUCTION 



The study of nonequilibrium phase transition univer- 
sality classes is one of the most fundamental tasks of 
statistical physics [1-6] . Such phenomena may appear in 
models of population, epidemics, catalysis, cooperative 
transport [1] and enzyme biology [7] for example. For a 
long time only the universality class of the directed per- 
colation (DP) [8] has been known, but later it turned out 
that different classes may appear in other models (for a 
review see [3]). The most prominent example occurs in 
one-dimensional systems exhibiting Z 2 symmetrical ab- 
sorbing states [9-13] and parity conserving (PC) - ex- 
plicit or underlying - branching and annihilating random 
walk (BARWe) dynamics (A -> 3A, 2A -> 0) of the dual 
variables (A) [14-17]. For a review see [18]. This class 
is also called as directed Ising, DP2 or generalized voter 
model class. 

However understanding these phenomena involving a 
satisfying solution of the corresponding field theory is 
very rare. For BARWe models field theory failed to give 
quantitatively precise results in one dimension [17] be- 
cause systematic epsilon expansion brakes down due to 
a second critical dimension at d! c = 4/3 below d c = 2. 
Very recently an other field theory has been suggested 
and analyzed by numerical simulation of the Langevin 
equation for systems exhibiting Zi symmetric absorbing 
states [19]. Numerical approximations ranging from sim- 
ulations [13,15,18,20,21] to series expansions [22], clus- 
ter mean-field approximations [11,21] and empty interval 
method [23] established the values of critical exponents 
firmly. 

Recently Zhong et al. [23] constructed a special parity 
conserving reaction-diffusion model (hereafter we call it 
ZAM model). By studying it up to pair approximations 
they claimed, that cluster mean-field studies fail to re- 
produce the phase transition of the BARWe model if one 
considers clusters which are not large enough. Therefore 
one may arrive to the conclusion that cluster mean-field 
and especially the site mean-field solution break down 
generally in case of PC class transitions. However in 
that study the branching attempt probability (a) was 
fixed to a finite value, so one should not expect to see 
the mean-field transition, which is known to occur at 
zero branching rate. In the present study we show that 



the cluster mean-field approximations can describe the 
the mean-field transition qualitatively well even for small 
cluster sizes if we do not exclude the neighborhood of the 
zero branching rate from the parameter space. These ap- 
proximations are performed on one-dimensional clusters 
therefore for large cluster sizes N one expects to see a 
convergence towards the PC class transition at a c > 0. 
We apply coherent anomaly extrapolations for the se- 
quence of cluster mean-field results and give estimates 
for some exponents of the critical behavior in one dimen- 
sion. 

The interplay of diffusion and fluctuation has already 
been shown in many reaction-diffusion models (see for 
example [21,27,29,30]). Although the diffusion is unable 
to change the universal behavior it can affect the loca- 
tion of the transition and even more it can change the 
stability of a fixed point in case of competing reactions 
[27,29,30]). To investigate the possible role of diffusion 
we also extend the parameter space by modifying the 
diffusion rate. 



II. THE CLUSTER MEAN-FIELD METHOD 

The generalized (cluster) mean-field (GMF) method is 
an extension of the usual mean-field calculation by set- 
ting up master equations for n-point configuration prob- 
abilities of site variables Sj € {A, 0} 



dP n ({sj}) 
dt 



f{Pn{{Si\)) 



(1) 



where the function / depends on the transition rules of 
{Pm] "block probabilities" at time t. 

At the level of A-point approximation the correlations 
are neglected for n > N, that is, P n (si, s n ) is ex- 
pressed by using the Bayesian extension process [38-40] , 



Pn (^1 j Sn) 



HjJo N Pn(si+3,- 



. s N+j) 
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(2) 



In principle, 2 N — 1 parameters are required to define 
the probability of all the Appoint configurations. This 
number, however, is drastically reduced by the following 
conditions : In the stationary state the particle distribu- 
tion is assumed to be symmetric with respect to trans- 
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lation and reflection. Furthermore, the block probability 
consistency results in: 

P„(si, s n ) — P„ + i(si, s n , s„+i) , 

s n +l 

P n (si, S n ) = 2J Pri+l( s 0, Si, S n ) . 

so 

Here we apply GMF for one-dimensional, site re- 
stricted lattice versions of BARWe. Taking into account 
spatial the symmetries in case of the N — 10 GMF ap- 
proximation one has to find the solution of equations of 
528 independent variables. This has been achieved with 
the help of MATHEMATICA software. We required 20 
digit accuracy in the results and arbitrary precision dur- 
ing the calculations. 

It is well known that such approximations predict the 
phase structure qualitatively well in one dimension, pro- 
vided N is large enough to take into account the rele- 
vant interaction terms. For example N > 1 is needed to 
take into account particle diffusion terms, while N > 2 
was found to be necessary in case of binary production 
processes involving pair induced reactions [27,28]. The 
GMF is an efficient phase diagram exploration method 
and although it is set up for the d = 1 lattice in previ- 
ous cases it provided qualitatively good phase diagram 
for higher dimensional, mean-field versions as well (see 
example [26,29,35,36]). 

III. THE ZAM MODEL 

In [23] Zhong et al. defined a special, one-dimensional, 
parity conserving lattice model in which each site is ei- 
ther empty or singly occupied. The state of the lattice is 
updated asynchronously by the following rules: an occu- 
pied site is chosen randomly and it is tried for diffusion, 
at rate T (probability T/(T + Q)), or branching, at rate 
fi (probability Q/(T + £1)); while the time is increased 
by 1/N, where N is the number of occupied (active) 
sites. In a diffusion step the particle jumps to its ran- 
domly chosen nearest neighbor sites. If the site is oc- 
cupied both particles are annihilated with probability r. 
On the other hand the jump is rejected with probabil- 
ity 1 — r. The branching process involves the creation of 
two new particles around the neighborhood. If cither, or 
both neighboring sites is previously occupied the target 
site(s) become empty with probability r. Otherwise the 
lattice remains unaltered with probability 1 — r. 

The site mean-field equation for the concentration c 
(probability of site occupancy) is 

4c = -2rc 2 + 2c(l - c) 2 - 2rc 3 , (3) 
dt 

where r = fi = 1 was taken, hence the branching proba- 
bility is fixed to 1/2. Similarly the N — 2 pair mean- 
field approximation was solved and in both cases the 



system always evolves to an active phase with finite con- 
centration c s (see Fig.l). These results were compared 
with those of the r < 1 rate model in one dimension, 
for which PC class continuous phase transition is known 
at r c > [15]. The conclusion was drawn that cluster 
mean-field approximations of low orders fail to repro- 
duce the phase diagram, the convergence is very slow 
and the calculations are complicated. However it is also 
known that in the mean-filed approximation - which is 
valid for d > d c = 2 dimensions - the phase transition 
occurs at zero branching rate [17], so it is not surprising 
that by fixing fl = 1 one does not see a phase transition 
in the corresponding mean- field approximations. On the 
other hand Zhong et al. found a phase transition for 
< r c < 1 by the parity interval method - which also 
involves a mean-field like approximation - set up for this 
one-dimensional model [23]. The authors acknowledged 
that it had been possible to see the phase transition in 
related PC class models for larger (N > 2) cluster sizes 
[11]. 

To clarify this we extended the cluster mean-field 
method for higher orders and for other versions of this 
BARWe model exhibiting PC class transition. 
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FIG. 1. (Color online) GMF results for the steady state 
concentration in the ZAM model for N = 1, 2, ..7 clusters 
(top to bottom). 



We determined the steady state solutions of Eq.(l) tak- 
ing into account Eqs.(2,3) and calculated the correspond- 
ing steady state densities c s (r) (see Fig.l). Indeed higher 
levels of approximations (N > 4) of the ZAM model re- 
sult in phase transition with r c (N) < 1 converging to- 
wards the simulations (r c = 0.470(5)). 
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IV. GMF RESULTS FOR THE ZAMB MODEL 



c s oc \r - r c \ 



(7) 



As we saw in the preceding section cluster mean- 
field approximations of the ZAM model give qualita- 
tively good phase diagram for N > 4 and the steady 
state solutions converge to the MC simulations. How- 
ever going much further with the GMF study of the ZAM 
model is time consuming especially because the numeri- 
cal root finding of Eq.(l) gets computationally demand- 
ing for large number of variables. On the other hand 
by modifying the ZAM model slightly in such a way 
that zero branching rate is allowed with the restriction 
a = 2ft/ (T + Q) = 1 - r (ZAMB) one immediately finds 
the expected mean-field transition at a c = for N > 1. 
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FIG. 2. (Color online) Steady state concentration GMF re- 
sults for the ZAMB model for N = 1,2,3,4,5,6,7,8,10,12 
clusters (top to bottom). 



and (3' = 2 correction to scaling exponent defined as 

c s = a\r-r c f + b\r-r c f . (8) 

The steady-state solution has been determined for 
N = 1,2,3,4,5,6,7,8,10,12 (see Fig.2). As one can 
see even the pair approximation gives r c (2) = 1, but 
for N > 2 the transition point starts shifting towards the 
true transition point r c < 1 as expected in one dimen- 
sion. One can also observe a concave shape of the curves 
corresponding to the corrections to scaling with [3' = 2 
exponent. 
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FIG. 3. (Color online)Steady state concentration GMF 
results for the ZAMB model for D — 0.2 diffusion and 
N = 1,2, 3, 4, 5, 6, 7, 8, 10, 12 clusters (top to bottom) 
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The site approximation equation for the concentration 

is 

— c = 2c (l + c 2 (r-l) 2 -r-c (2- rj) , (4) 
at V / 

exhibiting the steady state solution 

2-r- yW8- llr + 4r 2 

Cs = 27^T? • (5) 

for r < 1. For r = 1 Eq.(4) simplifies to 

|c = -2c* (6) 

resulting in c oc l/t particle density decay in agreement 
with the mean-field expectations. The supercritical be- 
havior in the active phase can be characterized by f3 = 1 
leading order singularity for r c < 1 



A. The effect of diffusion 

In previous papers [11,25-27,29,30]) it was shown that 
the diffusion strength can cause relevant effects on the 
phase diagram of reaction-diffusion models when it com- 
petes with the reactions. This is reflected in the cluster 
mean-field approximations in such a way that stronger 
diffusion "washes out" fluctuations and causes a transi- 
tion, which is more site mean-field like, while higher 
takes into account more fluctuations, hence opposes the 
effects of the diffusion. Here we investigated the effect 
of diffusion by lowering the hopping probability of the 
ZAMB model to D = 2T/(T + O) = 0.2. The steady 
state results (Fig. 3) show that the concentration curves 
arrive with higher slopes to the c s = axis than in case 
of D = 1. This permits us to obtain a(N) more pre- 
cisely since the relative error of the amplitudes is smaller 
and the quadratic correction to scaling is weaker. As the 
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consequence the numerical root finding is not so much af- 
fected by basin of attraction of the absorbing state fixed 
point solution. Note that in the pair approximation the 
linear amplitudes a(2) are zero. 

V. SIMULATION RESULTS 

By applying fitting with the expected scaling form 

\r c (N) - r c \ v ± oc l/N (9) 

to the GMF data one can determine the location of the 
transition and u± simultaneously. For D — 1 this gives 
r c = 0.402, while for D = 0, 2 r c = 0.65. However to ob- 
tain better critical value estimates we used more precise 
r c values in the fitting procedure, which can be deduced 
from simulations. 

1.1 i 1 1 1 — i 
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FIG. 4. (Color online) Particle density decay of the ZAMB 
model times the expected critical power-law (i ' 285 ). Differ- 
ent curves correspond to r = 0.565, 0.562, 0.56 and D — 0.2, 
r = 0.408, 0.409, 0.41 and D = 1 (from top to bottom). 



b) If i — Sj-|-i = A the two nn particles are removed 
with probability r. 

c) If s,_i =/= s i+1 the nn are swapped with probability r. 
The time (t) is updated by l/np again. The simulations 
were followed up to t = 10 7 (MCS) or until n P = (ab- 
sorbing state). The concentration of particles c s (t) times 
the expected density decay power-law 1/c oc £ - 285 [18,20] 
is plotted on Fig. 4. One can read-off r c = 0.409(1) for 
D = 1 and r c = 0.562(1) for D = 0.2. For D = 1 this 
agrees well with the extrapolation results of (Eq.(9)), but 
for D = 0.2 the deviation is not negligible. 

VI. COHERENT ANOMALY EXTRAPOLATIONS 

According to scaling theory the location of the critical 
point for sizes N (r c (N)) scales (in the large N limit) as 
Eq.(9). Precise extrapolation can be obtained by apply- 
ing the critical transition point values of simulation in the 
mentioned scaling form. The r c (N) for the N-th level of 
approximations was determined by quadratic fitting for 
c s (N) < 0.002. Figure 5 shows r c (N) as the function 
of l/N. The fit of form (9) yields l/v ± = 0.54(1) for 
D = 1 and l/i/± = 0.53(1) for D = 0.2 (see Fig.5). The 
value v±_ = 1.85(3), agrees well with the value from the 
literature v± = 1.84(6) [20] for this class. 




The simulations were performed on one-dimensional 
lattices of sizes L = 10 with periodic boundary con- 
ditions. The runs were started from half filled lattices 
with randomly distributed particles. One elementary 
Monte Carlo step (MCS) consists of the following pro- 
cesses. A particle and a direction are selected randomly. 
If the nearest neighbour (nn) in the selected direction was 
empty the particle moves to it (with probability D). If it 
was filled, both particles are removed (with probability 
r). The time - measured by MCS — is updated by l/np, 
where np is the total particle number at time t. To per- 
form the branching another particle is selected randomly 
(with probability 1 — r). Depending on the status of the 
two nn (sj_i,Sj_)_i) the following process may occur, 
a) If Si-x = Si + i = 0, two new particles are created. 
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FIG. 5. (Color online) CAM extrapolation results for the 
critical point r c (N) (2 < N < 12). The insert shows |o(iV)| 
(3 < N < 12). Stars correspond to D = 1, boxes to D = 0.2. 
Numerical errors are smaller than the symbol sizes. 



We applied the coherent anomaly method (CAM) [31] 
for the A^-cluster GMF results to extrapolate to the 
N — > oo behavior. This method has been proven to be 
successful for obtaining critical exponents of nonequilib- 
rium absorbing phase transitions [11,21,32-35,37]. Ear- 
lier the GMF+CAM method was already used for an- 
alyzing the PC class transition of the NEKIM [11,21] 
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up to N < 6 cluster sizes. That study arrived to the 
rough estimate (3—1, which agrees with the the series 
expansion result — 1.00(5) [41] and the simulation re- 
sults marginally (3 — 0.95(2) [18]. Now we apply this 
method for cluster sizes up to N < 12 and improve the 
GMF+CAM results [11,21] for a model PC exhibiting 
PC class transition. According to CAM the amplitudes 
a(N) of the cluster mean-field singularities 



c s (N) = \a(N)\\r c (N) 
scale in such a way 

|a(A0| a |r c (JV) - r c f~^ 



/3a 



(10) 



(11) 



that the exponent of the true singular behavior (Eq.(7)) 
can be estimated. The a(N) amplitudes were determined 
by linear fitting to the local slopes of the c s (N) data in 
the neighborhood of r c (N). The amplitudes are shown 
in Table I and in the insert of Fig. 5. The fitting using 
form Eq.(ll) for N > 2 data (since a(2) = 0) results in 
(3 = 0.92(5) for D = 1 and (3 = 0.96(2) for D = 0.2. 
The CAM for D = 1 results in bigger numerical error 
than for D = 0.2, because for D = 1 the quadratic cor- 
rections to scaling are stronger and the c s (r) numeri- 
cal solutions are affected by the attractive basin of the 
neighboring absorbing state fixed point. These values 
agree well with the simulation results (3 = 0.95(2) [18], 
(3 = 0.94(6) [20] as well as with that of the parity interval 
method = 0.92(2) [23]. 
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3 


0.9575(2) 


0.4124(4) 


0.8997(2) 


0.094(5) 


4 


0.9177(1) 


0.4155(3) 


0.8321(1) 


0.100(8) 


5 


0.8802(4) 


0.4167(9) 


0.7805(1) 


0.103(7) 


6 


0.8485(1) 


0.4192(9) 


0.7423(9) 


0.104(7) 


7 


0.8235(1) 


0.4254(3) 


0.714(7) 


0.100(12) 


8 


0.8044(1) 


0.429(3) 


0.691(13) 


0.103(9) 


10 


0.7787(1) 


0.427(2) 


0.6608(1) 


0.100(6) 


12 


0.7633(2) 


0.416(2) 


0.6394(1) 


0.103(6) 


r c 


0.562(1) 




0.409(1) 






1.88(3) 






1.85(3) 







0.96(2) 




0.92(5) 



TABLE I. Summary of results for the ZAMB model 



VII. CONCLUSIONS 



We showed that even low order GMF approximations 
can describe the phase transition of an even-offspringed 
BARW model correctly if one allows appropriate param- 
eterization. We applied GMF approximations for a one 
dimensional site restricted lattice model. By allowing 
to reach the zero branching rate in the ZAM model we 
showed that the site mean-field solution is in agreement 
with that of the field theory for this class [17]. Note that 
this kind of analysis resulted in similar steady-state solu- 
tions in case of an other PC class model [11,21], although 
there the branching rate can't be read-off explicitly. 

The GMF approximations were determined up to N — 
12 and convergence towards the simulation results were 
shown. Using scaling and CAM theory we obtained 
v± = 1.85(3) and = 0.96(2) critical value estimates 
matching the best precision available in the literature for 
the PC universality class. 
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